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Abstract. The small sizes of low mass stars in principle provide an opportunity to find Earth-like 
planets and "super Earths" in habitable zones via transits. Large area synoptic surveys like Pan- 
STARRS and LSST will observe large numbers of low mass stars, albeit with widely spaced (sparse) 
time sampling relative to the planets' periods and transit durations. We present simple analytical 
equations that can be used to estimate the feasibility of a survey by setting upper limits to the number 
of transiting planets that will be detected. We use Monte Carlo simulations to find upper limits for 
the number of transiting planets that may be discovered in the Pan-STARRS Medium Deep and 3n 
surveys. Our search for transiting planets and M-dwarf eclipsing binaries in the SDSS-II supernova 
data is used to illustrate the problems (and successes) in using sparsely sampled surveys. 
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PLANETS ORBITING LOW MASS STARS 

Planets transiting low mass (M < 0.6 Mq) stars are important for many reasons. To 
~10% accuracy, the radii of stars with M < 1 Mq are Rstar/R® ~ Mstar/MQ. Conse- 
quently, it is easier for ground based telescopes to detect planets of all sizes, especially 
Neptune-sized ice-giant planets, around M-dwarfs. Table [U illustrates this point by show- 
ing the ratio of geometrical transit depths expressed as magnitudes of Jupiter, Neptune, 
and the Earth transiting stars of different masses. The radii of the stars were calculated 
using Rstar/R® = 0.927 x Mrta,./M(^-|- 0.037 ±0.02, which is a linear fit to the masses 
and radii of 38 low mass stars (see UlSSBBl^ J])- 

With a limit of ~2 mmag light-curve precision, transiting Neptune-sized planets are 



TABLE 1. 


Radii, Spectral type, and geometric transit depth in 


magnitudes of transits 


for planets as 


a function of the mass of the parent star 








Mstar/M0 


Rstar/R0 


Spectral Type 


Mi*/M,. 


Jupiter 


Neptune 


Earth 


1.00 


1.00 


G2 




0.011 


0.001 


0.0001 


0.65 


0.64 


MO 


8.4/7.9 


0.028 


0.003 


0.0002 


0.45 


0.45 


M2 


8.9/8.3 


0.057 


0.005 


0.0004 


0.15 


0.18 


M4 


10.7/9.9 


0.447 


0.036 


0.0029 


0.092 


0.12 


M6 


13.2/11.9 


1.299 


0.076 


0.0061 


0.079 


0.11 


M8 


14.9/13.2 


2.118 


0.094 


0.0075 



very difficult to detect via transits around main sequence G-dwarf host stars, even in 
the best ground-based photometric data. However, a Neptune-sized planet orbiting a 0.4 
M0 M-dwarf has a detectable transit with a depth of 0.036 mag, and becomes more 
detectable as the host mass decreases. Observing these transits around M-dwarfs will 
yield the radii of Neptune-like planets and rocky planets. These radii, combined with ra- 
dial velocity measurements, determine the planets' densities, and thus provide important 
information about their mean composition. Reflected and emitted light measurements, 
in addition to transmission spectra, constrain the physics and composition of planetary 
atmospheres. 

ESTIMATING THE NUMBER OF PLANETS TRANSITING LOW 
MASS STARS IN A GIVEN SURVEY 

We will show why large area surveys are needed to find appreciable numbers of planets 
transiting low mass stars by beginning with some simple geometrical considerations. 
Defining the inclination of the planet's orbit as the angle between the line of sight and 
the orbital plane, the critical angle for a transit to occur is 

sine, = {R, + Rp)/R (1) 

where R^ and Rp are the radii of the star and planet respectively, and R is the radius 
of the circular orbit. If there are Np planet bearing stars with randomly oriented orbits, 
the number of stars with planets that can transit (as opposed to the number that will be 
observed), is 

No=Np{Rs + Rp)/R, (2) 

showing the strong selection for short period planets. The number of planets Nj transiting 
at any given instant is Np times the fractional area of a circle with radius Rs + Rp 
projected onto a sphere of radius R, and is given by the equation 

iV/ = A^pX(l/2)(l-cos0,). (3) 

Equation [3] is easily understood by thinking of the projection of Rs + Rp onto a sphere 
with radius R as a "polar cap", with the line of sight along the axis of the pole. 

We can now estimate the number Nt of planets that can possibly transit during a con- 
tinuous observation time that is less than half the planet's period by first considering 
an observation time ^1/4 that lasts for 1/4 of a period. If we suppose that all of the tran- 
siting planets have zero inclination, then half of the planets A^o can reach the polar cap 
and transit, and half of this number will have an orbit toward the polar cap, and in fact 
transit. Consequently, 

Nt^{Np/4)(R,, + Rp)/R. (4) 

Because not all transiting planets have zero inclination. Equation |4] is an upper limit to 
the number of transits that will be observed, and becomes a progressively better estimate 
as the critical angle 0c becomes smaller (i.e. as the period increases for a fixed stellar 
mass). Taking into account the fraction of planets that project onto the polar cap, an 



approximate upper limit for the number of planets that are transiting or will transit during 
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where the angle Qt is given by 

Qt = Oc + iTttobs/P^M. (6) 

We next estimate the number of stars N{DJ,b) that will be in a survey area to 
distance D corresponding to a limiting magnitude m/,m. Although the density of stars is 
a function of Galactic longitude and latitude, M-dwarfs are faint enough that, depending 
on miiin, only a few hundred parsecs will be sampled. Consequently, for a first order 
model, we can assume an exponential density as a function of Galactic latitude. The 
number of stars is then 

fD _£sini 

N{D,b)=No{mum2)da r^e dr (7) 

Jo 

where No{mi,m2) is the density of stars of a particular type in the galactic plane, and 
ro is the scale height perpendicular to the plane. To calculate NQ{mi,m2) we assume a 
mass function of the form 

^ = 0.0187m""stars pc"^ (8) 
dM 

where the normalization and coefficient a = 1.3 are taken from Reid et al. dsj]. The 
equation integrates to give 

0187 

A/o(mi,m2) = — [m\^" - m\^") stars pc"^. (9) 

We now use the preceding equations to estimate how many M-dwarfs of a given type 
will be in a particular survey. As an example, consider stars with masses between 0. 1375 
and 0. 175 Mq, corresponding to a mean spectral type of M4. Using Equation^ the space 
density of M4 stars is ~ 0.008 stars pc^-^; if the 5a limiting magnitude is z* = 21.0, the 
limiting distance is ~1700 pc, and there are ~1300 M4 stars/sq-deg at = 15 deg. 
However, we show in the next section that detection of Neptune-sized objects with a 
SNR of 5 in a single observation requires decreasing the 5c7 limiting magnitude by 3.6 
mags to z* = 11 A, dropping the number of M4 stars to ~300 per sq-deg. If we assume 
1 hour of continuous observing and that all of the planets are Neptune-like in a ~2 day 
orbit, the fraction of planets that can transit (Equation O is 0.05, and the fraction that 
will transit (Equations [5] and (6]) is 6.2 x 10^"^. If we make the optimistic assumption 
that 10% of all M4 stars have short-period planets, the probability that a transit occurs 
during a 1 hour observation is 6.2 x 10^^ per planet bearing star. We must survey nearly 
55 sq-deg for an hour to detect just one transiting planet! 

Surveying at low Galactic latitudes increases the number of M4 stars by an order of 
magnitude, at the price of decreasing the photometric accuracy because of crowding. 



and greatly increasing the number of false positives from blends of single stars with the 
target variable stars. Calculations show that observing at infrared wavelengths where the 
M dwarfs are much brighter helps very little from the ground because the gains from 
increased luminosities are completely offset by the increases in sky background. 

We look at the prospects of finding planets transiting low mass stars in sparsely 
sampled, wide area surveys below. 

THE PAN-STARRS-1 SURVEY 

The prototype telescope for the Panoramic Survey Telescope and Rapid Response Sys- 
tem (Pan- STARRSfi PSl, is presently being commissioned at the Haleakala Observa- 
tories in Hawaii. Three PSl programs are relevant for finding transiting planets. The 
first, Pan-Planets [2J, consists of 3 campaigns of z-band bright time imaging of several 
adjacent fields in the Galactic bulge, the Hyades, and in Praesepe. The Galactic bulge 
fields will include ~480,000 dwarf stars to a limiting magnitude of i = 17.0; ~100 short 
period transiting Jupiter-sized planets are expected to be found. The host stars will be 
bright enough for the planets to be confirmed by radial velocity observations. Unfortu- 
nately, the small number of M-dwarfs in these fields (~ 15,000; [@]) brighter than the 
limiting magnitude will lead to very few detections of Jupiter and Neptune-sized planets 
around these stars. 

The second of these programs, the PSl Medium Deep Survey (MDS), will survey up 
to 7 fields (one field is 7 sq-deg) every usable night, covering a total of 10 fields in a 
year. Each field will be observed in g and r, i, z, and y over the course of four nights. 
Assuming 65% of the nights are clear, the cycle of five filters in four nights will be 
repeated 35 times per year for three years. The 5c7 limiting magnitudes for an z-band 
exposure on a single night is given in Tabled The PSl goal for absolute photometry 
precision is 0.01 mag. To detect the 0.036 mag depth (0.033 in flux) of a Neptune 
transiting an M4 star in a single observation, the photometric noise in the star must 
be < 0.007 mags. Assuming that differential photometry will meet this requirement, the 
magnitude cutoff for the survey must be decreased by 3.6 magnitudes. The average of the 
estimated number of M-dwarfs with masses < 0.175 M© (~M3.5) in the MDS z-band 
and z-band images is ~ 24, 000. To first order the average of the numbers in the r and 
y filters will be comparable. The majority of the stars will have spectral types between 
M3.5 and M5. Taking M = O.1375M0 (~M4.5) as representative of these stars and 
assuming all of them have a Neptune -sized planet with a period drawn from a uniform 
random distribution of periods between 1 and 3 days, a Monte Carlo simulation of three 
years of observations gives the transits listed in Tabled 

Even if we assume that 10% of all late M-dwarfs have a short period Neptune-sized 
planet, the MDS numbers are relatively small, given that only a fraction of planets that 
do transit during an observation will be detected in the face of random noise, red noise, 
blending with eclipsing binaries or variable stars, and other effects. The good news is 
that there should be ~ 10 transits per short period transiting planet in the MDS survey. 



See: |http : //panstarrs . if a . hawaii . edu/publlc/| 



TABLE 2. Estimated upper limits to the number of detectable planets transiting low mass stars in the 
Pan-STARRS MDS and 3>n surveys. Columns 1-6 describe the assumed cadence for the MDS and 37r 
surveys, the 5(7 limiting magnitudes for a single exposure, and the 140(7 limits for a 5(7 detection of 
a Neptune-sized planet transiting an M4 dwarf. The last three columns give the number of stars in the 
respective survey areas to the 140(7 limit, the number of planets that will transit if 10% of these stars 
have a Neptune-sized planet in an orbit with a period drawn from a uniform random distribution between 
1 and 3 days, and the average number of transits per star 



Survey 


p., Num. Obs./ 
FielcWear 


Obs. 
Sep. 
Nights 


Survey 
Limit* 
(5(T)~ 


Survey 
Limit 

(140(7) 


Num. of 
Stars 


Trans. 
PlanetsJ. 


Avg. Num. 
Trans/Star 


MDS 


r,i,z,y - 140 


1 


25.4 


21.8 


-24,000 


< 120 


10.0 




i 2 


30 


23.35 


19.7 


9.7 X 10'' 


< 4,000 


1.04 


3k 


z 2 


30 


23.35 


18.7 


8.3 X 10*^ 


< 3,400 


1.04 



* The 3(7 limiting magnitudes in the 3;r survey have been increased by 0.75 mag based on the assumption that the 
two images in a single night that are separated by 30 to 60 minutes are combined. 

^ This number is derived by assuming that 10% of the stars in the Monte Carlo simulation have a short period 
Neptune-sized planet 



which will increase their detectability. 

Table [U shows that transiting Jupiter-sized planets are far more detectable around the 
earlier and brighter M-dwarfs than Neptune-sized planets. Gaudi fio'] concludes that 
~40 such planets may be found by the MDS. This estimate, however, depends on the 
unknown frequency of such higher mass planets around M-dwarfs. 

The final PSl program that we consider, the 3Tt steradian survey, will observe the 
entire sky visible from Hawaii approximately once a week. Each field will be observed 4 
times a year in each of the five filters for 3 years (60 epochs in all filters). Each night that 
a field is observed in a particular filter, the field will be observed twice with a separation 
of 30-60 minutes to identify moving solar system objects. We (arbitrarily) assume that 
the observations two nights per year for each filter per field will be separated by 30 days, 
and the transits observed in the i and z images are for different stars. A Monte Carlo 
simulation spanning three years yields the numbers given in Tabled Again, assuming 
that 10% of all low mass M-dwarfs have short period Neptune-sized planets, the 3k 
survey could yield ~ 300 such stars with 2 possible transits each. Detecting planets 
which transit only twice during widely separated short observation will be very difficult. 
Unless the cadence between the r,i,z,y filters is optimized for detecting transits in the 3ti 
survey, few planets transiting low mass stars will be detected. 

Another issue that must be addressed is stellar activity in fast-rotating, fully convec- 
tive low mass dwarfs. Two recent encouraging photometric campaigns on the young 
active MlVe star AU Mic demonstrated that the star's flare activity, its 4.85 day rotation, 
and variability due to two large star spots, did not prevent differential photometry accu- 
rate to a few mmag flT]. A modest photometric study of late M-dwarfs could be used 
to assess the limits on transit detection set by stellar activity. The next step is to make 
more realistic simulations that include noise sources, specific detection algorithms, and 
search for optimal survey cadences. 



THE SDSS-II SUPERNOVA SURVEY 



The SDSS-II Supernova Survey 1121], in its three years of operation, monitored an ~300 
sq-deg area on the Celestial Equator known as SDSS Stripe 82 for three months per 
season, with a rough cadence of one observation every two nights. We used photometric 
catalog^ from the Survey to search for low mass M-dwarf binaries and transiting 
planets. We extracted objects that were classified as stars with clean photometry H and 
enforced an upper magnitude limit of z = 21.0 mag. We then matched detections of 
objects across all three seasons to build up an extensive light-curve catalog. 

From this catalog, we selected objects that had > 10 observations over the three sea- 
sons, leaving us with approximate^ 1.1 million point sources. We used ensemble pho- 
tometry techniques lfl3ll to generate differential magnitude light-curve^, thus decreas- 
ing the chance of incorrectly flagging an object as periodically variable if its light-curve 
was affected by the systematic artifacts caused by changing sky conditions and/or small 
changes in photometric zero points. We obtained an average photometric precision of 
about 30 mmag at / = 19.0 mag, which effectively rules out the possibility of finding 
transiting planets around low mass stars, but is comfortably within the detection regime 
of M-dwarf eclipsing binaries. 



We used SDSS r - i and i - z color cuts defined by West et al. UM to select ~690,000 



likely M-dwarfs from our sample of 1.1 million point sources limited to z = 21.0 mag. 



To find periodic variables, we employed the Stetson variability index [|15|1 : this removed 
from consideration a large number of faint objects that appeared to be variable because of 
photometric noise. We then inspected the differential magnitude light-curves of color- 
selected likely M-dwarfs that survived this selection process, looking for multi-band 
events that were eclipse-like. 

A rough estimate of the total number Ndet of eclipsing binary M-dwarf systems we 
expect to find is given by 

Ndet =N^xfMX fbin (P) X ftrans X f thresh ■ (10) 

A^* is the total number of stars in our sample, is the fraction of M-dwarfs seen in 
the sample, fhm{P) is the M-dwarf binary fraction as a function of period, f trans is the 
fraction of M-dwarf binary systems that can be observed by our survey, and f thresh is the 
fraction of M-dwarf eclipsing binary systems that will be detected given our survey's 
SNR and throughput limits. Using the numbers from our sample for A^* ~ 1.1 x 10^, 
/m — 0.60, and conservative estimates for /fo„(2.5 — 7 days) ~ 0.01, f trans — 0.05, and 
fihresh — 0.05, we estimate that we may be able to find up to 17 M-dwarf eclipsing 
binaries. 

To date, we have processed ~30% of our sample and have found 12 candidates. We 
present the light-curves of two such objects in Figure [IJ We are currently testing various 
period-finding techniques, including BLS [fl6ll . and AoV [17], to determine the likely 



periods of these systems. This will allow us to generate ephemerides for photometric 



Available at: |http : / /www . sdss . org/ dr sn 1 /DRSNl_data_re lease . html] 

|http : //cas . sdss . org/ dr5/ en/help/ docs/algorithm. asp| 

See M. Richmond's implementation: [http :// spiff . rit.e du/ ensemble/] 



MB5920, type: Ml, 41 observations 
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FIGURE 1. Differential photometry riz light-curves for two M-dwarf eclipsing binary candidates from 
the SDSS-II Supernova Survey: MB5920 (left panel) and MB15389 (right panel). The differential magni- 
tudes have been shifted on the y-axis for clarity. Both candidate names are internal designations. 



monitoring, and enable eventual radial velocity foUowup to confirm whether they are 
indeed eclipsing binary systems of low mass M-dwarfs. 
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